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ABSTRACT 

We present the results from the first large (> 100 source) 3.5 mm polarimetric survey of radio loud active 
galactic nuclei (AGN). This wavelength is favorable within the radio-mm range for measuring the intrinsic linearly 
polarized emission from AGN, since in general it is only marginally affected by Faraday rotation of the electric 
vector position angle, and depolarization. The /, Q, U, and V Stokes parameter observations were performed with 
the XPOL polarimeter at the IRAM 30 m Telescope on different observing epochs from July 2005 (when most 
of the measurements were made) to October 2009. Our sample consists of 145 flat-radio-spectrum AGN with 
declination > -30° (J2000.0) and flux density > 1 Jy at ~ 86 GHz, as measured at the IRAM 30 m Telescope 
from 1978 to 1994. This constraint on the radio spectrum causes our sample to be dominated by blazars, which 
allows us to conduct new statistical studies on this class of high-luminosity, relativistically-beamed emitters. We 
detect linear and circular polarization (above minimum 3cr levels of ~ 1.5 %, and ~0.3 %) for 76 %, and 6 % of the 
sample, respectively. We find a clear excess in degree of linear polarization detected at 86 GHz with regard to that 
at 15 GHz by a factor of ~ 2. Over our entire source sample, the luminosity of the jets is anti-correlated with the 
degree of linear polarization. Consistent with previous findings claiming larger Doppler factors for brighter 7-ray 
blazars, quasars listed in our sample, and in the Fermi Large Area Telescope Bright Source Catalog (LBAS), show 
larger luminosities than non-LBAS ones, but our data do not allow us to confirm the same for BL Lac objects. 
We do not find a clear relation between the linear polarization angle and the jet structural position angle for any 
source class in our sample. We interpret this as the consequence of a markedly non-axisymmetric character of 
the 3 mm emitting region in the jets. We find that intrinsic circular polarization is the most likely mechanism for 
generation of the circular polarization detected in our observations. Our new data can be used to estimate the 
3.5 mm AGN contribution to measurements of the linear polarization of the cosmic microwave background, such 
as those performed by the Planck satellite. 

Subject headings: Galaxies: active - galaxies: jets - BL Lacertae objects: general - quasars: general - 
polarization - surveys 



1. INTRODUCTION 

Radio loud active galactic nuclei (AGN) are known to pro- 
duce powerful pairs of highly collimated relativistic jets of 
magnetized plasma, which are able to extend far beyond 
the boundaries of the host galaxy. Their relativistic nature, 
their magnetic fields, and their non-isotropic geometry deter- 
mine their most rele vant observational properties: superlu - 
minal motions (e.g., iGomez et al. 1 12001 1 iJorstad et al. 1 12005b . 
Doppler boosted (decreased) emission of the jet pointing at a 
small (large) angle to the observer (e.g., Kadler etaTl 120041). 
rapid intrinsic emission variability (e.g. 
iFuhrmann et al.ll2008HMarscher et al. 



Agu do et al.1 12006b: 



and changes of the 



jet structure (e.g.. IJorstad et al. 20071 Agudo et alj|2007b . and 
intense polarized synchrotron and inverse-Co mpton emission 
along all the electromagnetic spectrum (e.g ., lOstorero et al.l 
120061 lAbdo eTaD2010):lMarscher et al.ll2010h . 

Polarimetric very long baseline interferometry (VLBI) ob- 



servations at mm wavelengths, combined with single dish ra- 
dio, mm, and optical linear polarimetric observing campaigns, 
allow one to (i) connect the location of the emitting regions 
at different observing bands and (ii) infer properties about the 
nature of the innermost moving and stationary knots of emis- 
sion in the jets and the magnetic field in s uch regions, which 
usually cannot be resolved by VLBI (e.g.. IJorstad et al.l 120071: 
iMarscher et al.ll2008ll20Toh . 

VLBI polarimetric surveys of radio lo ud AGN also provide 
relevant information about their jets; see iPollack et al.l J2003) 
for a survey on a sample of 177 sources in the Calteck-Jodrell 
Bank Flat-Spectrum survey observed with the Very Long Base - 
line Array (VLBA) at 5 GHz, and Lister & Homanl d2005l) : 
iHoman & Listed (120061) for the first results from the MOJAVE 
survey on 133 s o urces obser ved with the VLBA at 15 GHz. 
iLister & Homanl (120051) and IPollack et al.1 d2003l) found the 
source cores (i.e., the innermost visible jet region at a given ob- 
serving frequency with an instrument capable of resolving such 



1 Instituto de Astrofi'sica de Andaluci'a (CSIC). Apartado 3004, E-18080 Granada, Spain 

2 Institute for Astrophysical Research, Boston University, 725 Commonwealth Avenue, Boston, MA 02215; 

iagudo@bu.edu 

3 Institut de Radio Astronomie Millimetrique, 300 Rue de la Piscine, 38406 St. Martin d'Heres, France; 

thum@iram.fr 



4 Instituto de Radio Astronorm'a Milimetrica, Avenida Divina Pastora, 7, Local 20, E-18012 Granada, Spain; 

wiesemey@iram.es 

5 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel, 69, D-53121 Bonn, Germany; 

tkrichbaum@mpifr-bonn.mpg.de 



1 



2 



Agudo et al. 



a jet) to be weakly linearly polarized ( 5 %), but with signif- 
icantly larger fr actional polarizat i on for the jet regions down- 
stream. Indeed, iLister & Homanl d2005l) reported a general in- 
crease of linear polarization degree with increasing distance 
outward from the core in quasars and BL Lacertae (BL Lac) ob- 
jets, but with BL Lac jets more polarized in general. They also 
found the cores and jets of radio galaxies to be much weaker or 
not line arly polarized compared to those in quasars and BL Lac 
objects. iHoman & Listed ( 120061) found "strong" circular polar- 
ization (>0.3%, as defined by them), usually in the cores of the 
sources (but not always), within ~15% of their sample. They 
did not find significant correlations between the degree of circu- 
lar polarization in the core and other relevant source properties. 

No polarization survey of a large number of AGN has been 
performed at mm wavelengths thus far. It is known that AGN 
jets and their cores are affected by non-negligible Faraday ro- 
tatio n measures (RM) typically ranging from ~ 10 2 radm~ 2 
(e.g.. IZavala & Tavlorf l2003l l2004t iGabu zda et al.ll200U l2004t 



Asa da et al.l 120081) to ~ 10 4 radm" 2 (e.g. jZaval a & Tavlorl 



2002tlAttridge et alJ2005UG6mez et alJ2008h . These large RM 



values ~ 10 4 radm~ 2 have only been found for a few sources , 
typically in the jets of radio galaxies dZavala & Tavlorl [20041) . 
Since the observed electric vector linear polarization angle of 
a source Xobs = Xmt + RM A 2 (where RM A 2 is the amount of 
Faraday rotation and Xmt is the intrinsic polarization angle), a 
A =3.5 mm (86 GHz) survey would be much less affected by 
Faraday rotation than are cm wavelength observations. At this 
wavelength, |RM| ps 7000 rad m" 2 would be required to produce 
Faraday rotation by ps 5 °, the median linear polarization angle 
uncertainty in the measurements presented in this paper. In con- 
trast, such a value of |RM| translates into rotations of ps 160° 
and ps 1400° at 2 cm (15 GHz) and 6 cm (5 GHz), respectively. 
If AGN jet RM in the radio and mm emitting regions are similar 
in general - which is still to be confirmed observationally- the 
reduced effect of Faraday rotation in mm wavelength observa- 
tions also makes them much less sensitive to Faraday depolar- 
ization, hence reflecting the intrinsic linear polarization proper- 
ties of the sources with better fidelity than radio observations. 

The mm wavelength emission of radio loud AGN is typi- 
cally dominated by compact regions in their jets; either in their 
innermost core s of emission, or in their pc scale super lumi- 
nal knots (e.g., iJorstad et alj|2007l: iMarscher et al.ll2008b. In- 
deed, the results from the 86 GHz VLBI survey of Lee et al.l 
d2008h have shown that, in general, ~ 50 % of the total flux 
density of their imaged sources is concentrated in regions not 
further than ~ 1 .5 milli-arcseconds from the mm core. In con- 
trast, because of the longer lifetime of synch rotron radiation at 
cm wavelengths dRvbicki & Lightmanlll979l) . such emission is 
more spread out al ong the pc scale of the jets up to distances > 
5 milli-arcseconds dPollack et alf2 003: Lis ter & Homanll2005h . 
Moreover, owing to synchrotron opacity, the bright cm wave- 
length emission from compact cores in AGN lies in regions 
farther downstream in the jet than at mm wavele ngths (e.g., 
lLobanovll 19981: iMarscheifeOOfl IJorstad et alj|2007l) . 

In the absence of a large mm- VLBI polarimetric survey, the 
statistical analysis of a single-dish polarization-mm survey can 
yield relevant insights into the intrinsic polarization properties 
of the innermost compact regions of relativistic jets in AGN. 
Here we present the results from the first 3.5 mm polarimetric 
survey over a large (> 100 source) sample of radio loud AGN. 
Our results are also a valuable resource for surveys dedi- 



cated to studies of the polarization properties of the cosmic mi- 
crowave background (CMB), specially that conducted by the 
ESA's Planck mission (see the mission BlueboolQ). For such 
studies, radio loud AGN are considered the most important 
CMB-foreground contributors away from the Galactic plane at 
angular scales ^ 0.5° up to ~ 100 GHz observing frequency 
dTucci et all 12005b . Hence, a priori knowledge of their po- 
larization is important to account for their contribution to the 
polarization of the CMB d ata (e.g.. lRubifio-Martm et alj|2010h 
lLopez-Caniego et"al]|2009b . 

2. THE SAMPLE 

In Table Q] we present our sample of 145 observed sources 
(see also Fig. [TJ. This sample has been built essentially from 
the 138 radio loud AGN in the IRAM pointing source list. Such 
list comes from the sample of ~ 300 radio-to-mm flat-spectrum 
co mpact AGN v i sible f rom the IRAM 30 m Telescope selected 
by ISteppe etaL J 1988b that were contained in the 1 Jy cata- 
log dKiihr et aljl 19811) and/or in the catalog of positions , struc- 
tures, and polarizations of 404 compact radio sources bv lPerlevI 
(fl982h . 

From this sample. ISteppe et al.l d 199211 19931) and lReuter et al.l 
d 1997b selected the brightest 138 AGN sources measured by the 
IRAM 30 m AGN monitoring program with a typical flux den- 
sity S90 £ 1 Jy at ~ 90 GHz and a minimum at S90 ~ 0.5 Jy in 
the time span from 1978 to 1994. After removing the 4 sources 
with J2000.0 declination < -30° (too low to be observed at the 
30 m Telescope above 20° elevation), we complemented our 
sample with 1 1 new sources not contained in the IRAM point- 
ing source catalog but strong enough currently to be considered 
within the selection criteria for such a cat al og. 

Since the samp le by ISteppe et alj d 19921 1 1993b and 
iReuter et alj d 1997b was essentially defined to be flux density 
limited at 5go £ 1 Jy from 1978 to 1994, our sample was ex- 
pected to be complete in that sense. Most likely because of 
source variability, 54 (19) sources in our sample have measured 
flux densities 5s6 < 1 Jy {S^ <0.5Jy), hence not matching, at 
the observing epoch, the constraints imposed for the 1978 to 
1994 time span. 

To have a raw estimate of the number of radio loud AGN 
with ~ 1 Jy that we may have missed in our sample, we 
compared it with the 5-Year WMAP Point Source Catalog 
(WMAP5. IWright et al.ir2009h . From this catalog, we find 28 
sources with J2000.0 declination > -30°, and with flux density 
larger than 1 Jy at 90 GHz but not contained in our sample. In 
contrast, our sample contains 20 sources with 5s6 > 1 Jy that are 
not found in WMAP5. All but one of these are low Galactic 
latitude sources, which reflects the bias of the WMAP5 catalog 
that excludes sources near the Galactic plane. 

Assuming that a 1 Jy flux-limited complete sample of radio 
loud AGN above -30° dec. contains all sources with 586 > 1 Jy 
both in our sample (i.e. 91), and in WMAP5, we only miss 
24 % of eligible sources. Thus, although this does not allow 
us to claim completeness of our sample to the 1 Jy total flux 
density limit, the relatively small fraction of missing sources, 
together with the unprecedented large size of our 86 GHz po- 
larization sample, allows us to rely on the results derived from 
it as long as subsample selection does not lead to unacceptably 
small number statistics, which is not the case for any of our 
subsamples (see below). 

In a survey like ours, polarization non-detection is a relevant 
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bias that may be introduced by source variability when the total 
flux density of a considerable fraction of the sample decreases. 
However, we have checked that we are not very much affected 
by such bias. Among our 91 sources with S&(, >1 Jy, 77 % were 
detected in linear polarization, whereas for the 54 sources with 
$86 < 1 Jy> the fraction of linear polarization detections is simi- 
lar (74 %). 

Given the spectral criteria applied to select the sources for the 
IRAM pointing source list, our sample is dominated by radio- 
to-mm flat-spectrum compact AGN, i.e., by blazars. To be more 
specific, it contains 107 quasars, 26 BL Lac objects, and 6 radio 
galaxies, with 6 unclassified sources, i.e., not contained in the 
I Veron-Cettv & Veronl d2006l) catalog. The sampl e also contains 
32 so urces in the Large Area Telescope (LAT, lAtwood et al.l 
120091) 7-ray bright AGN source sample (LB AS). The LBAS list 
contains sources detected by LAT onboard the Fermi Gamma- 
ray Space Telescope durin g its first 3 months of operation (i.e., 
those in Tables 1 and 2 of lAbdo et ai1l2009l) . These 32 LBAS 
sources include 19 quasars, 12 BL Lacs, and, 1 radio galaxy . 
Note that out of the 117 LBAS sources in lAbdo et ail (120091) . 
only 97 are visible from the 30 m Telescope with declination 

> -30°. Among these 97, 32 (33%) are in our 3.5 mm sam- 
ple. The source redshift in the sample ranges from z = 0.004 
to z = 3.408, with mean and median z = 0.941 and z = 0.815, 
respectively. 

Hence, our sample is adequate for studies of the mm polari- 
metric properties of quasar and BL Lac blazars, and of their 
relation with their 7-ray properties. Both the MOJAVE and 
our samples were selected primarily based on the flatness of 
source radio spectra. Hence, both of these are dominated by the 
sa me kind of sources, i.e., blazars. Among the 133 sources in 
the lLister & Homanl d2005l) MOJAVE sample, 94 (71 %) are in 
our sample as well. Both samples, and the results derived from 
them, can hence be directly compared, as we do throughout this 
paper. 

3. OBSERVATIONS AND DATA REDUCTION 

The observations w ere performed with the XPOL polarime- 
ter dThum et al.ll2008l) on the IRAM 30 m Telescope, by mak- 
ing use of the Observatory's orthogonal linearly polarized A100 
and B100 heterodyne receivers tuned at 86 GHz (3.5 mm). The 
main observing block was performed on July 2005. To com- 
plete the list of observed sources, a small number of measure- 
ments had to be performed in a time block in September 2005 
or at different epochs between 2006 and 2009 within different 
observing programs; see Table|2] 

All observations were performed under the standard XPO L 
set-up and calibration scheme discussed in lThum et al. I (l2008h . 

Every XPOL polarization measurement consisted of a se- 
ries of wobbler switching on-offs with typical integration times 

> 4min., depending on the source's total flux density. Every 
polarization integration was preceded by a cross-scan pointing 
of the telescope and an amplitude, phase, and decor relation loss 
polar ization-specific calibration measurement (see iThum et ail 
120081 

To estimate the remaining instrumental polarization to be 
subtracted from the data (a posteriori), we made measurements 
of strong, compact and unpolarized sources at almost every ob- 
serving epoch. For 86 GHz observations at the 30 m Telescope, 
the HPBW of the beam is ~28", hence an object with an ap- 
parent size of ( ^ 4") can be safely considered as compact for 
our purposes. As a rule, we used Mars, Uranus, and Neptune 



as instrumental polarization calibrators. When not available, 
other planets with up to ~12" at the time of the observation 
(e.g., Venus), the compact HII regions W3 0H, K3-50A ,and 
NGC 7538, and the planetary nebula NGC 7027, were observed 
for cross-check. 

The output of every polarization observation consists of a 
wide band (640 MHz) spectral measurement from both the 
A 100 and B100 receivers, together with the real and imagi- 
nary part of their cross-correlation. These four observables are 
needed to recover t he 4 Stokes param eters for each measure- 
ment as described in lThum et al.l (l2008h . After the observations, 
the polarization calibration measurements were used to cali- 
brate the amplitudes of the A 100 and B 100 measurements, their 
relative phase, and the decorrelation losses. Such calibrations 
were applied within the MIRA and CLASS software packages 
inside GILDAS. The instrumental polarization for every ob- 
serving epoch (\Qi\ < 2%, \Ui\ £ 0.5 %, and \Vi\ < 0.5 %) was 
then estimated and removed from the data. The Cj y /K = 6.4 Jy/K 
calibration factor for 86 GHz o bservations at the 30 m Tele- 
scope (e.g.. lAgudo et al. I l2006bl) was also applied to the total 
flux density measurements (S^)- At this stage, the uncertainties 
in every specific measurement only account for statistical un- 
certainties on the average of data along the spectral bandwidth. 
The final errors in the total flux density measurements were 
computed by adding quadratically a 5 % systematic fact or com- 
ing from the uncertainties in Cj y /K dAgudo et al.ll2006bl) . In this 
stage, the remaining polarimetric non-systematic errors, which 
determine our polarimetric precision, were estimated from the 
dispersion in the Q, U, and V Stokes parameters from our main 
unpolarized calibrators, i.e. Mars and Uranus. These disper- 
sion estimates (Agi =0.5 %, At/; =0.3 %, and AVj =0. 1 %) were 
also added in quadrature to provide the final uncertainty for ev- 
ery polarization measurement. The latter translate into final 
polarization-uncertainty medians of Awl ~0.53%, Ax ~5°, 
and Amc «0.2 % for the linear polarization degree (mjj, the 
linear polarization electric vector position angle (x), and the cir- 
cular polarization degree (mc), respectively. The uncertainties 
of our best measurements are dominated by AQi, At/;, and AV;, 
and hence their final uncertainties are Atol «0.5 %, Ax ~1.5 °, 
and Am c «0.1 %. 

4. RESULTS AND DISCUSSION 

In Table [2] we present the observing epoch and integration 
time for every source in our sample, as well as the observational 
results, expressed in terms of 5g6, tol, X> an d m c- A 3a upper 
limit in both m^ and mc is given whenever the measurement did 
not exceed such 3a value. Note that whereas linear polarization 
was detected from most sources in our sample (76 %), circular 
polarization could be detected only for a small fraction of them 
(6 %). 

Statistical analysis and discussion of the relevant aspects re- 
garding these data are presented in the following sub-sections. 

4.1. Total Flux Density 

4.1.1. Total Luminosity 

Fig. |2] shows the 86GHz luminosity (L = 47raf£S 86 (l + z)~\ 
where di is the luminosity distance for a H„ = 71 km s" 1 Mpc -1 , 
ft m = 0.27 and S7a = 0.73 cosmology, used hereafter) as a func- 
tion of redshift for our sample, with most — all but 19 sources 
— above the Sg(, = 0.5 Jy threshold in the observer's frame. 
Only 6 sources in our sample displayed S&(, > 5 Jy. As expected 
by the dependence of L on cosmological distance, quasars, the 
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most distant AGN class, show the largest luminosities (with me- 
dian L Q = 3.1 x 10 27 W/Hz), followed by BL Lac objects (L B = 
2.7 x 10 26 W/Hz), and radio galaxies (L G = 3.8 x 10 24 W/Hz), 
the closest AGN (see also Fig. 0). 

The LBAS sources in our sample show an L distribution 
which is apparently different from that of non-LBAS sources 
(Fig. [3j. The median and peak luminosities of LBAS sources 
(Llbas = 2.4 x 10 27 W/Hz, L%£ s within [3,10] x 10 27 W/Hz) 
are larger than those for non-LBAS ones (L non -iBAs = 1.9 x 
10 27 W/Hz, L^-lbas within [1,3] x 10 27 W/Hz). However, the 
Kolmogorov-Smirnov (K-S) test does not give a sufficiently 
high confidence level (i.e., 70.7%) to conclude that both dis- 
tributions are selected from different parent distribution^ 

Comparison of the luminosity of 7-ray classes of quasars 
and BL Lacs separately (Fig. [U shows that LBAS quasars 
in our sample are apparently more luminous (LLBASquasars = 
7.0 x 10 27 W/Hz, L^ B a A k SqHamrs within [3,10] x 10 27 W/Hz) 
than non-LBAS quasars (Lmn-LBASquasars = 2.7 x 10 27 W/Hz, 
L P nZ-LBAS„ within [1,3] x 10 27 W/Hz), whereas BL Lacs 
show similar median luminosities and luminosity distributions 
independently of their 7-ray properties. Indeed, the K-S test 
shows that the distributions of LBAS and non-LBAS quasars 
are significantly different (within a 99.0 % confidence level of 
both being drawn from different parent distributions), whereas 
LBAS and non-LBAS BL Lacs are not (confidence level only 
13.8%). 

Note that, for the computation of the confidence level of re- 
jection of the null hypothesis in the K-S test, the effective num- 
ber of data points (/V^fi must be > 4 to rely on the K-S test 
dPress et al.lll992l) . All K-S tests presented in this paper, in- 
cluding the one for LBAS and non-LBAS BL Lacs mentioned 
a bove (with N e = 5.7 ), fulfill this requirement. 

iLister et al.l d2009l) found that among the 26 MOJAVE AGN 
contained in the LBAS sample, quasars have, on average, 
faster superluminal features than non-LBAS quasars. They 
also report evidence of faster superluminal proper motions in 
7-ray variable LBAS AGN with regard to non variable ones. 
iKovalev et alj d2009f) report a correlation of the 7-ray photon 
flux of 77 LBAS sources with contemporaneous VLBA or sin- 
gle dish flux density at 15 GHz. They also show that LBAS 
sources observed in their 15 GHz VLBA sample display, in 
general, larger compact radio flux densities, larger brightness 
temperatures of their parsec-scale cores, and more active radio 
states than non-LBAS sources. These results are in good agree- 
ment with previous findings based on EGRET data suggest- 
ing that 7-ra y bright blazars posses l a rger Doppler factors than 
weak ones (Jorstad et all 2001alb1 : lLahteenmaki & Valtaoial 
120031: iKellermann et al.l l2004t) . However, ILister et al 
also suggest that the relativistic Doppler factor is not the sole 
parameter controlling the 7-ray properties of blazars. 

Our results are consistent with previous claims suggesting 
that 7-ray bright blazars have larger Doppler factors than weak 
7-ray blazars, but we can only confirm the larger luminosity of 
7-ray bright quasars against 7-weak quasars, and no significant 
difference for BL Lacs. 

We note, though, that our 86 GHz observations were not per- 
formed contemporaneously with the Fermi-LAT 7-ray obser- 



vations bv lAbdo et alj (|2j)09j), and the luminosity distributions 
of LBAS and non-LBAS sources presented here (specially for 
BL Lacs) may appear broader by 3 mm source variability, thus 
hiding actual trends (see § |4.4t . 



4. 1 .2. 75 GHz to 86 GHz Spectral Index 

The spectral index between two observing frequencies 
(o„,,„ 2 = ^ g(SujS V2 )/log{v\/v2)) in the radio-mm spectral 
range provides information about the synchrotron opacity be- 
tween such observing frequencies {vi,vz) affecting the radia- 
tion from an emitting jet region. For the 15 GHz and 86 GHz 
emission in radio loud AGN, the emitting regions are not co- 
spatial in general (see Section[T]i. Hence, for most sources, the 
0:15,86 determination is biased towards larger (more optically 
thick) spectral indices. This is because the innermost 86 GHz 
jet emitting regions are usually opaque to 15 GHz radiation. 

With this in mind, we performed an 015,86 study of our main 
subsamples only to study the relative differences among them. 
For this, we used only those sources in our sample with avail- 
able 15 GHz total flux densities from integrated intensities of 
MOJAVE VLBA imagefl To try to avoid biasing the computa- 
tion of 015,86 owing to source variability, we selected, for each 
source, the MOJAVE observation closest in time to our 86 GHz 
measurement. 
The 

015,86 distributions for each one of the source samples 
considered in this paper are shown in Fig. [5] Even when af- 
fected by the above mentioned bias, the spectral indices for the 
whole source sample are distributed towards flat and optically 
thin spectral indices (with 015,86 median 015,86 = -0.23, as ex- 
pected from the definition of the sample and the high observing 
frequency), with a small fraction of sources (18%) showing 
flat to optically thick spectral index. Quasars and non-LBAS 
sources posses 015,86 distributions similar to those of the entire 
source sample {af^ = -0.27, <5™ n 8 - LBAS = -0.26). However, 
both BL Lacs and LBAS sources clearly distribute spectral in- 
dices more uniformly over an almost symmetric range of 015.86 
in [-0.7,0.6]. This points out that BL Lac and LBAS sources 
tend to show considerably flatter mm spectra (partially optically 
thick) than quasars and non-LBAS sources, respectively. The 
significance of this result is guaranteed by the K-S test, which 
points out that the 015,86 distributions for quasars and BL Lacs 
in Fig. [5] are different (at 99.4% confidence level), as well as 
those for LBAS and non-LBAS subsamples (within a 96.9 % 
confidence level). We note that quasars in the LBAS subsample 
are dominated by sources with spectral indices above Oj""^ 1 , 
whereas LBAS BL Lacs have a more spread 015,86 distribution 
along the ranges of the entire BL Lac subsample. 

In principle, the spectral differences shown by radio loud 
quasars and BL Lacs in our analysis would be expected if jets in 
BL Lacs would be preferentially oriented towards smaller view- 
ing angles (with regard to the line of sight). In this case, the in- 
tegrated radiation coming from differ ent spectral compo nents is 
expected to flatten t he spectrum (e.g., Cotton et aljl980h. How- 
ever, recent work dHovattaetalJl2009r Pus hkarev et al. I l2009h 
has concluded that flat spectrum radio quasars have signifi- 
cantly smaller viewing angles than BL Lac objects, which is 
expected to produce the reverse spectral behavior from that ob- 
served by us. In this case, we can attribute more confidently 
the average spectral differences shown by quasars and BL Lacs 

7 As usual, only confidence levels > 95 .0 % will be considered sufficiently high to claim that two distributions are significantly different throughout this paper. 

8 N e = N\Nz/(,N\ +N2), where JV,- is the number of points on each sample, as defined by Press et al. 1 1992) 

9 |http : / /www . physics ■ purdue ■ edu/MOJAVE| 
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in Fig. [5] to the larger cosmological redshifts of quasars (see 
Fig.0, which shift their spectra to lower frequencies in the ob- 
server's frame. This preferentially reveals the optically thin part 
of the synchrotron spectrum of quasars with regard to BL Lacs. 
Also, we cannot rule out differences in the intrinsic spectral 
properties of the emitting particle populations between these 
two source classes. 

The means of the viewing a ngles of LBAS and non -LBAS 
sources contained in Table 1 of iPushkarev et a n d2009h are not 
statistically significant according to their Student's T-test. This 
does not allow us to attribute the spectral differences between 
these two subsamples to preferential alignment of 7-ray bright 
blazar jets with the line of sight, although this possibility can- 
not be ruled out based on only this information. Intrinsic differ- 
ences in the properties of the synchrotron-emitting particle pop- 
ulations of these two subsamples, as well as source redshifts, 
can, in principle, be relevant as well. 

4.2. Linear Polarization 

Fractional linear polarization at 86 GHz (mi) was detected 
for 110 sources, 76 % of our sample. When comparing the me- 
dian values of tol (wl) for the major optical classes, we find that 
BL Lac objects, with = 4.4%, are more strongly polarized 
than quasars, with tol = 3. 1 %. To avoid biassing the hil values, 
3a,„ L upper limits were considered for the cases when tol could 
not be detected. Otherwise, if the non-detections are not ac- 
counted for when calculating the medians, tol would be overes- 
timated by ~ 1-1.5 The difference between the quasar and 
the BL Lac distributions is confirmed by the Gehan's gen- 
eralized Wilcoxon (GGW) test at a 95.2 % confidence level. To 
perform this test, which takes into account both detections and 
uppe r limits, we used the ASURV 1 .2 survival analysis package 
(see lLavallev et alJI 19921 and references therein). 

This result (i.e., quasars less polarized at 86 GHz, in gen- 
eral, than BL Lacs) does not seem consistent with quasars hav- 
ing a significantly thinner synchrotron spectrum than BL Lacs 
between 15 and 86 GHz. This is expected to reduce for 
BL Lacs with regard to quasars owing to stronger synchrotron 
self-absorption in BL Lacs. Hence, the radio-to-mm spectral 
properties of quasars and BL Lacs cannot explain their over- 
all linear polarization properties. In contrast, an explanation 
comes from recent evidence that the vie wing angle of jets in 
quasars is smaller tha n that in BL Lacs dHovatta et al.l 120091: 
IPushkarev et al.ll2009l) . If either the magnetic field is not ho- 
mogeneously distributed along the jet or the jet has prominent 
non-axysimmetry, lower polarization degree is expected from 
sources better oriented to the line of sight (i.e., quasars) owing 
to cancellation of orthogonal polarization components. 

Interestingly, the otl distribution of the entire source sample 
is double peaked (Fig. [6]), with the first peak at otl w 2.5 %, 
and the second one at otl « 4 %. At lower polarization de- 
gr ees, a similar bimodal d istribution was previously observed 
bv lLister & Homanl d2005l) for the cores of quasars and the in- 
tegrated polarization degree of EGRET-detected blazars. Fig. [6] 
shows that this double peak in our data might come from the 
quasar subsample, which dominates the overall sample. 

To test whether this dichotomy is produced by the pres- 
ence of radio loud quasars with high optical polarization (HPQ, 
> 3 %) and radio quasars with low optical polarization (LPQ, 
< 3 %) in our sample, we have also analyzed these two sub- 

10 Note that assigning different values to the non-detected does not have an 
as is the case for all considered subsamples. 



samples separ ately (see Fig.[6]l. We consid ered HPQ as those 
sources in the lVeron-Cettv & Veronl d2006l) catalog with such 
a label, and LPQ as the remainder of sources in that catalog. 
Simple inspection of Fig. [6] apparently shows that if there is a 
physical meaning on the two-peaked quasar distribution, it can- 
not be because of HPQ producing the higher peak and the 
LPQ the one at lower m^. Indeed, there is no statistically sig- 
nificant difference between the HPQ and LPQ distributions 
that could confirm such a hypothesis (the GGW test gives only 
a 21.4% confidence). 

A similar double peak is seen in the ;«l distribution for the 
LBAS source sample (Fig.|6]l, whereas the second peak is not so 
evident in the non-LBAS source hil distribution. This might be 
interpreted as suggesting that this bimodal distribution comes 
from 7-ray bright quasars. However, according to our GGW 
test, there is no significant difference between the otl distri- 
butions of LBAS and non-LBAS sources (36.9 % confidence 
only). 

Thus our analysis remains inconclusive about the origin of 
the apparent polarization degree dichotomy. 

4.2.1. 86 GHz to 15 GHz Fractional Linear Polarization Ratio 

The fraction of linear polarization detections in our sam- 
ple is <~ 76%. Considering that we adopted a 3a criterion 
for detection (with median 3<t,„ l ^ 1.6 % over the whole sam- 
ple), this means that, on average, ^76 % of our sources display 
wl ~ 1 .6 % at 86 GHz. However, only 60 % of sources both in 
the MOJAVE and in our sample show 15 GHz fractional linear 
polarization ^ 1.6%. This points to a general trend of blazars 
to increase their linear polarization degree w ith frequency, as 
previo usly suggested bv lAgudo et al.l d2006bl) and lJorstad et al.l 
Indeed, those sources for which linear polarization was 
detected both in o ur survey and in the inte grated emission of the 
MOJAVE images dLister & H oman 2005) follow a clear general 
trend to show significantly larger fractional linear polarization 
at 86 GHz than at 15 GHz by a median factor w 2 (Fig.|7). The 
same factor w 2 is obtained if our tol data is compared with 
the linear po larization degree of the core at 15 GHz (m™\ e 5 ) as 
measured bv lLister & Homanl d2005l) (Fig. O. 

Note also that there is a prominent 22 % (24 %) fraction of 
the entire source sample with larger tol/otl.is (mi/m^ 5 ) ra- 
tio than 4. This tail of large 86 GHz linear polarization ex- 
cess seems to be present in all subsamples of optical and 7-ray 
classes, although BL Lacs do not show it clearly, perhaps be- 
cause of the lower number of sources in that subsample. 

This result may be explained by a combination of two phe- 
nomena, that the 86 GHz emission in blazars comes from a re- 
gion with greater degree of order of the magnetic field than the 
one at 15 GHz, and that the 15 GHz emission from blazars is af- 
fected by considerably greater Faraday depolarization relative 
to the 86 GHz emission. 

If the former hypothesis is true, it would imply that the mag- 
netic field order in the inner regions of the relativistic jets is 
larger than in outer regions. For this statement, we assume that 
the bulk of the 86 GHz emission comes on average from inner 
jet regions where 15 GHz emission is strongly affected by syn- 
chrotron op acity, hence locating the 15 GHz core farther down- 
stream (e.g. Jjorstad et alj|2007l) . 

Despite the apparent differences between the source subsam- 
ples presented in each of Fig.|7]and[8](see also above), our K-S 

on the resulting median, provided that all estimates of oil are smaller than m^, 
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test does not allow us to claim any significant differences be- 
tween the m^/m^s or m^/m^W distributions either for quasars 
versus BL Lacs (confidence level only 84.5%, and 53.9%, 
respectively), or for LBAS against non-LBAS sources (confi- 
dence level only 85.1 %, and 59.5 %, respectively). 

4.2.2. Total Luminosity vs. Linear Polarization 

Fig. [9] shows the 86 GHz luminosity vs. 86 GHz fractional 
linear polarization for sources with known redshift in the entire 
source sample and in the four major optical and 7-ray classes 
considered here. Our correlation analysis for data containing 
upper limits — perf ormed with the ASURV 1.2 package by 
iLavallev et all d 19921) — shows significant correlation between 
L and otl for the whole source sample. This is supported by the 
results of the Cox, Kendall's r, and Spearman's p tests. These 
indicate that L and wil are correlated at 98.8%, 99.2%, and 
98.7 % confidence level, respectively. Significant correlation 
is also found for the case of the non-LBAS subsample — at 

98.4 %, 97.0 %, and 96.7 % confidence levels, respectively. 
However, for the quasar subsample — with confidence levels 

of 87.7 %, 93.3 %, and 91.1 % for the Cox, the Kendall's r, and 
the Spearman's p tests, respectively — , no formally significant 
correlation is found. Perhaps because of their decreased sizes, 
the correlation tests over the BL Lac and the LBAS yield con- 
fidence levels of correlation of these subsamples that are even 
much lower. 

Under the standard assumption that, in the absence of strong 
Faraday depolarization (which we assume to be the case for our 
3 . 5 mm observations, see SectionQ]), higher degree of linear po- 
larization reflects greater magnetic field order, we have found 
that the magnetic field order — in our entire source sample — 
increases with decreasing mm luminosity. 

A reasonable hypothesis to explain this phenomenon would 
simply involve orientation and relativistic effects. In principle, 
those sources whose jets are better oriented to the line of sight 
are expected to display larger luminosities (because of their 
larger Doppler factors) and also lower linear polarization de- 
grees (because of cancellation of orthogonal polarization com- 
ponents along the line of sight). To test this h ypothesis, we 
used t he "variability Doppler factors" given by iHovatta et alj 
( 20091 to compute the beaming corrected luminosities of the 
70 sources in their and our sample. The correlation analysis 
of all source types in Fig. [9] points out that, although the re- 
duced number of sources in every subsample tends to reduce 
the correlation between L un b e amed and otl, the decreasing trend 
of Lunbeamed with increasing tol still seems to be present. De- 
spite the reduction of data points, such anti-correlation is now 
statistically significant for the LBAS subsample (with 99.1 %, 

99.5 %, and 98.4% confidence level, from the Cox, Kendall's 
t, and Spearman's p test, respectively), which do not allow us 
to support the beaming scenario as a reliable explanation for the 
L and tol anti-correlation. 

This may have a plausible alternative interpretation that, al- 
though still qualitative and speculative, may be a valid proposal 
to investigate further. The AGN jet regions responsible for the 
3 mm emission are thought to be dominated by plasma dynam- 
ics, in contrast to the innermost Poynting flux-dominated jet re- 
gion, wher e the collimation and acceleration is assumed to t ake 
place (e.g.. Uorstad et alj|2007t iMarscher et alj|2008l 1201 Ol) . It 
is known that increasing jet kinetic energies with regard to their 
surrounding low speed wind favor the generation of fluid in- 
stabilities, leading to turbulence, in the region separating the 



inner-fast jet and the oute r-slow wind (e.g jMizuno et al]|2007h 
iMeliani & Keppensl2009l for different numerical setups for this 
scenario). This picture fits into our observing results. A faster 
jet would be more luminous, but would also produce more tur- 
bulence in the shear layer where the external wind penetrates 
the jet, hence reducing the magnetic field order and the linear 
polarization degree. 

4.2.3. Linear Polarization Angle versus Jet Position Angle 

In Fig.[l0]we show the distribution of misalignment of linear- 
polarization electric-vector position-angle (%, given in Table [2| 
with the jet structural position angle (</>j et , given in Table [TJ, 
i.e., \x~ $jet|> for the entire source sample, and the subsamples 
of quasars, BL Lacs, LBAS, and non-LBAS sources. 

Different angular resolutions achieved at different VLBI 
observing frequencies probe different jet regions that may 
show, in some part icular cases, large <f>p t differences (see e.g., 
lAgudo et all 120071 for the case of NRAO 150 with A0 jet « 
100°). Opacity effects of the inner jet regions at low frequen- 
cies can also impede a reliable determination of |x - 0jet| if </>jet 
is measured at v << 86 GHz. To try to avoid these biases as 
much as possibl e, we first searche d for <j}- jet values in the 86 GHz 
VLBI Survey bv lLee et al l d2008l) . When 86 GHz VLBI images 
were n ot available, the 15 GH z data from either the MOJAVE 
survey dLister & Homanl2005[ preferentially ) or those from the 
2 cm VLBA Survey dKellermann et al.ll2004l) were used. Other- 
wise, a deeper search was done for every source from references 
1 to 10 on Table [TJ The higher frequency result thus found was 
then adopted. 

Fig. [10] shows a weak trend in the entire source sample of 
alignment of \ close to 0°-30° to 0j et , such that \ lies almost 
parallel to the jet axis. A similar pattern is found in the quasar 
and non-LBAS subsamples. In contrast, the values of |x-</>j e t| 
in the BL Lac and LBAS subsamples tend to lie within the 
ranges 10°-30° or 70°-80°. However, our K-S tests indicate 
that there is no significant difference between the quasar and 
the BL Lac distributions (only 18.7 % confidence to come from 
the same parent distribution), and between the LBAS and non- 
LBAS distributions (only 63.9 % confidence). 

We therefore find no clear trend in any of the source sam- 
ples considered in this work for x to be aligned either parallel 
or perpendicular to 0j et . Actually, even for the entire source 
sample, which shows an apparent preference of aligned \ par- 
allel to </>jet, there is only a small excess of sources (~25 %) 
distributed at |x~ 0jet| < 30°. Similar results are obtained from 
the quasar and LBAS distributions. Even this ~ 25 % excess 
must be interpreted with care. When only highly polarized 
source states (m^ > 3 %) are considered - which seems to re- 
veal more clearly the y to 4>\ sl relation in the case of the quasar 
3C 454.3 (Uorstad etal] 120101) -, such excess is decreased to 
~ 18%. Moreover, Fig. QT] shows that, if the |x — $jet| distri- 
butions are generated from a smaller number of sources (only 
those brighter than 1 .5 Jy at 86 GHz), such an excess disappears 
completely. All this shows that if there is a preference to align \ 
parallel to the jets of the sources shown in Fig.[l0] such a pref- 
erence is intrinsically weak at 86 GHz and it is also partially 
hidden by random changes in the |%— </>j e t| distributions. 

Theoretically, the electric vector position angle (x) of the 
linearly polarized emission from relativistic axisymmetric jets 
should be o bserved either parallel or perpendicula r to the 
jet axis (e.g, iLvutikov etal] 120051: ICawthornd 120061) . How- 
ever, the results from several observational attempts to con- 
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firm this b i-modality do not show a robust agreeme nt among 
each other. iGabuzda. Pushkarev. & Cawthornd ( 120001) . through 
5 GHz VLBI observations, reported a tendency of the core and 
jets in a set of 25 BL Lac objects to show \ either parallel (pref- 
erentially) or^ergejidjcuhujo the jet direction. However, on the 
one hand, iPollack et al.l d2003l) found a strong tendency for the 
cores of their observed quasars to possess values of \ that are 
perpendicular to the jet axis, whereas no correlation was found 
for the jets in quasa rs, nor for the jets or core s in BL Lac objects. 
On the other hand, iLister & Homanl d2005l) found the cores of 
sources in the MOJAVE sample to show \ values aligned pref- 
erentially parallel to the jet axis, which was found to be a par- 
ticularly strong tendency for the case of BL Lac objects. Such 
a preference was not found in the jets of their sample, although 
BL Lacs did show a better tendency to have \ aligned with the 
jets. The apparent lack of consensus from these studies might 
be the result of th e well known frequenc y dependent Faraday 
rotation of x (e.g. JZavala & Tavlorll2004l) . 

Our 86 GHz polarimetric observations show that there is not 
a clear trend in our data for |x~0jet| to be distributed either near 
~ 0° or near ~ 90°. This diffe rence between our results and 
those bv lLister & Homanl d2005l) cannot be attributed to the use 
of (/>j et measured at 86 GHz in our analysis. Indeed, for the 40 
sources in our sample having measurements of both 0j e t.>43GHz 
and </>jet,i5GHzi both magnitudes are in general good agreement 
(with Spearman's correlation rank p = 0.91 at 99.9% confi- 
dence, and mean of </>jet,>43GHz _ </ , jet,i5GHz ~ 9°). Also, when 
</>jet,i5GHz is used instead of </>j e t,>43GHz for these 40 sources, the 
distribution in Fig.[l0]does not show noticeable changes. 

The median (maximum) uncertainty in the determination of 
X for those sources in our sample detected in linear polarization 
is S\ ~ 4° (<5x max ~ 9°). Thus, these uncertainties cannot ex- 
plain the lack of a strong trend in our |%— </>j et | distributions for 
X to be aligned either parallel or perpendicular to <j>j et . 

The large amplitude and rapid linear polarization variability 



at mm (e.g ., Jorstad et alj l2007h and radio wavelengths (e.g., 



lAller et alf I2003L and references therein), for which blazars 
are well known, is certainly an unavoidable effect both in our 
|x — <f>jet\ study and in those by previous authors (see § 14.41 ). 
However, we do not find strong arguments supporting a sce- 
nario that hides a clear \\— (j>j et \ trend in our data when allow- 
ing for it in previous radio surveys. Considerably larger x vari- 
ability amplitude, or a larger probability to find sources in high 
variability states at mm wavelengths could certainly explain this 
scenario, but there is still no strong observational support for it, 
as far as the authors know. 

There is though a likely expla nation for the dif f erence s be- 
tween our re s ults a nd those by ILister & Homanl d2005l) and 
IPollack et alJ (l2003h . Our data, which are essentially free of 
Faraday rotation (see Section [T), suffer from prominent syn- 
chrotron opacity well upstream the 15 GHz and 5 GHz cores, 
hence reflecting the properties of inner jet regions than the sec- 
tions of the jet that radiate predominantly at radio frequencies. 
Apart from other biases, differences on the |% — </>j et | distribu- 
tions a mm and radio wavelengths could be reflecting the prop- 
erties of the regions from where most of the emission at such 
wavelengths comes from. 

There is an increasing number of AGN jets where large 
mm-x rotations (>> 90°) with typical time scales of days to 
months are d etected either in the core s or in bright moving 
features (e.g., iD'Arcangelo et al.l 120071: lLarionov et alj l2008t 

^Ihttp : / /www .rssd.esa.int /index. php?proiect=Planck| 



iJorstad et al.l l2010h . Successful explanation of these phe- 
nomena requires non-axisimmetric jet dynamics such as emit- 
ting features that propagate down the jet along helical paths 
dMarscher et al.ll2008ll2010l) . Gradients of axisymmetric Fara- 
day rotation screens along the path of movin g jet features 
may also pro duce large x rotations at 7 mm (e.g jGomez et al.l 
200C| 1200 ll) . but much larger than typical RM values (>> 
10 4 radm" 2 ) are required to produce rotations at 3 mm to ex- 
plain Figs. [IO]and[IT] 

4.2.4. AGN as Source of Contamination of the 86 GHz Linear 
Polarization CMB 

Fig. |6] show that linear polarization of flat spectrum radio 
loud AGN is larger than previously expected from centimeter 
wavelength surveys (e .g., |Pollack etaT1l2003l : iRicci et aLll2004t 
ILister & Horn an 2005). This is clearly indicated by the excess 
by a factor of ~ 2 of 86 GHz linear polarization relative to 15 
GHz linear polarization (Fig. [7]), which points out that simple 
extrapolation from measurements at cm wavelengths system- 
atically underestimates the CMB -foreground polarization con- 
tribution of AGN at 3 mm wavelength. The median fractional 
linear polarization of detected sources is ~ 4 % with a con- 
siderable population (27 %) showing linear polarization > 5 % 
(Fig.|6|). Fig.[T2lshows that, at 3 mm wavelengths, linear polar- 
ization in radio loud AGN may reach significant values relative 
to the CMB contribution for high sensitivity measurements. In- 
deed, 62 % (20 %) of our linear polarization detected sources 
show |2| or \U\ > 50mJy (> 100 mJy). 

Hence, our observing results in Table [2] are relevant, and 
can be used, to estimate the 3 mm AGN contribution to lin- 
ear polarization CMB measurements, such as those performed 
by the Planck satellitePl A new epoch of observations, over 
an improved flux density limited, complete sample, is already 
planned for mid 2010. This will provide contemporaneous 
measurements with those of Planck, as well as the possibility to 
study the effects of linear polarization variability on the CMB 
radio loud AGN foregrounds. 

4.3. Circular Polarization 

The distribution of \mc \ for all major samples considered in 
this paper is presented in Fig. [T3] For comparison, we show 
separately the distributions of \mc\ detections (> 3a), \mc \ re- 
sults > 2(T, and all |mc| measurements (independently of their 
significance). 

The most evident result with regard to circular polarization is 
the low level of detection in our observations (~ 6 %, 8 sources 
out of 145), which prevents us from conducting a statistical 
study of circular polarization. This low detection rate is not 
surprising given the known low degree of circular polarization 
of blazars at radio w avelengths (typically ; $ 0.5 % at 2 cm as 
observed with VLBI. IHoman & Liste rll200"6h and the minimum 
mc detection level in our observations (3er mc ^ 0.3 %). 

Our circular polarization dete ction rate is lower than half of 
that of IHoman & Lister! d2006l) (~ 15 %), which may be ex- 
plained by the difference in angular resolution between their 
15 GHz VLBA observations and ours. IHoman & Lister! d2006l) 
showed that most circular polarization is confined to the core of 
the sources. If this also applies to the 86 GHz emitting regions, 
our single dish observations (to which all such emitting regions 
contribute) must reflect lower mc and a lower detection rate for 



our similar sensitivity to circular polarization (er„ 



0.1%). 
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Among our 8 sources with detected circular polarization, 
only 3C 84 (0316+413 ) and 3C 454.3 ( 2251+ 158) were also 
detected at £ 3a by iHoman & Listerl ((2006). They found 
significant circular polarization in the core and in 3 down- 
stream jet emission regions in 3C 84, and in one jet region in 
3C 454.3. The amoun t and sign of mc measured in 2003 by 
IHoman & Listerl d2006l) in the core of 3C 84 is not consistent 
with our single dish me asurement in 2005. Note though that the 
VLBI measurement by IHoman & Listed d2006l) . who detected 
strong circular polarization with opposite sign than in the core 
within a nearby jet region, may not be comparable with our sin- 
gle dish measurement. In contrast, both mc and its sign in the 
core of 3C 454.3 as measured by Homan & Lister's, and our 
measurements are consistent, despite the extreme mm flare dis- 
played by this source duri ng our observations (e.g. jRaiteri et al.l 
2008| iJorstad et al.ll2010T) . These discrepancies, as well as the 
fact that our circula r polarization detections do not agree in gen- 
eral with those by IHoman & Listed (120061) . is consistent with 
the k nown circular pol arization variability in radio loud AGN 
(e.g.. lAller et alj|2003l) . which m ay differ considerab ly at dif- 
ferent observing frequencies (e.g.. lHoman et al.ll2009h . 

Circular polarization in AGN jets may be generated either 
by intrinsically circularly polarized synchrotron radiation, or 
by Faraday conversion o f line ar into circular polar ization - see 
IWardle & Homanl d2003l) and lHoman elaTI d2009l) for detailed 
reviews of the different magnetic field configurations and con- 
ditions suitable for both mechanisms. Both mechanisms have a 
dependence with the emitting wavelength (A), that conspires to 
reduce mc with decreasing A, hence providing an additional ex- 
planation for the low level of mc detection in our observations. 
However, whereas for the intrinsic mechanism mc mt oc A 1 / 2 , 
for Faraday conversion the dependence is m uch more stronger, 
mc F oc A 5 under the assumptions adopted bv lWardle & Honiara 
J20031 Hence, whereas rac 111 ' at 3 mm is only ~ 2 times lower 
than at 2 cm, mc ¥ is reduced by a factor ~ 6, 000. Additionally, 
some of the Faraday conversion scenarios invoked in the litera- 
ture dWardle & Homanl 120031: iRuszkowski & Begelmanl 120021: 



iBeckert & Falckd l2002h ~l Homan et al. 1 120091) involve Faraday 
rotation (internal to the jet), that has a A\ F oc A 2 depen- 
dence. All of these considerations render the Faraday conver- 
sion mechanism less suitable to explain the circular polarization 
detected at 3 mm. Otherwise it would be difficult to explain 
that the amount of circular polarization detected in our obser- 
vations, 0.3% ^ \mc\ 0.7%, is similar to that detected in the 
MOJAVE survey at 2 cm dHoman & Listedl2006l) . 

4.4. Influence of Source Variability 

It is well known that radio loud AGN display in the mm 
range (as in other spectral ranges) abrupt total flux density vari- 
ability by up to ma xima of ~ 1 order of magnitude and time 
scales from months ( Jorstad et al. 2005 : Terasranta et al.|[2005t 



iFuhrmann etafl 12008^ to days dAgudo et a lj |2006bl) . Such 
variability is usually connected with the ejection of superlu- 
minal emission fea t ures from the innermost jet regions (e.g ., 
IJorstad et all 120051: iKadler et al.l 120081: iPerucho et all l2008h . 
The Doppler factor of the source affects its variability time scale 
in the observer's frame through a directly proportional relation, 
whereas the cosmological redshift slows down such time scale 
b y a factor ( 1 +z)~ l ■ 

IJorstad et al.l d2005l) made a study of the (7, 3, and 1) mil- 
limeter polarization variability of 15 radio loud AGN moni- 
tored during ~ 3yr. They identify different classes of vari- 



able sources, with the most variable one displaying excursions 
of up to ~ 1 order of magnitude in linear polarization de- 
gree, and > 90° in linear polari zation angle in tim e scales of 
months or even weeks (see also IJorstad et alj|2010l for a case 
of mm and optical polarimetric correlation in 3C 454.3). At 
7 mm, shorter tol and \ variability time scales of less than a 
we ek - but for smaller variability amplitudes - are reported 
bv ID' Arcangelo et al.l (l2007l 120091) . for the case of the blazars 
PKS 0420-014, and OJ 287. 

However, not all radio loud AGN vary as much as in these 
extreme cases. There seems to be a gradient of variability am- 
plitude and time scales, which is likely related to the intrinsic 
properties of every source. Anyhow, source variability is an 
unavoidable effect that influences our results, as well as those 
from any other survey at which AGN display significant vari- 
ability. Such influence on the selection of our sample and its 
completeness was discussed in §[T] 

For the large number of objects in different variability states, 
our distributions of S&c, m^, and \ stre expected to broaden 
their shapes or to decrease their correlation between them and 
with other variables, hence conspiring to hidden possible sig- 
nificant results from them. However, for such a large sample 
like ours, we do not expect that variability can fake the sta- 
tistical results by shifting the distributions towards a particular 
trend. From that point of view, we are confident on the results 
in previous sections, which significance and extension to other 
subsamples could perhaps have been more ample without the 
influence of variability. This is extensive to the studies of 015.86 
and otl/otl.is in § |4.1.2| and |4.2.3l respectively. The extra vari- 
ability at 15 GHz that we could not account for introduces extra 
widening of the distributions against significant statistical re- 
sults. Even in that case, we were able to obtain some significant 
results. 

5. CONCLUSIONS 

We have performed the first large 3.5 mm polarization sur- 
vey of radio loud AGN over a sample of 145 bright and flat 
radio spectrum sources, dominated by blazars, observed with 
the IRAM 30 m Telescope. At such wavelengths, Faraday ro- 
tation and depolarization are expected to have a weak effect on 
the observed linear polarization, which thus reflects the intrinsic 
linear polarization properties of the sources with high fidelity. 

We detect linear and circular polarization above 3er levels for 
76 % and 6 % of the sample, respectively. We have shown that 
the fractional linear polarization at 86 GHz clearly exceeds that 
at 15 GHz by a factor of ~ 2 for all classes of sources consid- 
ered here. This implies both a larger degree of magnetic field 
order in the region where the bulk of the 86 GHz emission is 
produced, and considerable Faraday depolarization at 15 GHz, 
although we can not quantify the relative magnitude of these 
effects with the present data. 

Consistent with previous work suggesting that the class of 
7-ray bright blazars possess, in general, larger Doppler factors 
than weaker 7-ray bright blazars, we have found that LBAS 
quasars from our sample are more luminous than non-LBAS 
quasars. However, our data do not allow us to claim differences 
in the luminosity distributions of LBAS and non-LBAS BL Lac 
objects. 

Our entire source sample shows a trend of lower total flux 
luminosity for larger degrees of linear polarization, indicating 
that the level of magnetic field order in the innermost jet emit- 
ting regions is anti-correlated with jet luminosity. 
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Unlike other authors, who observed at radio frequencies 
(hence affected by Faraday rotation), we do not find a clear rela- 
tion between the linear polarization angle and the jet structural 
position angle of BL Lacs. This is also true for the remain- 
der of the source subsamples considered in this work, despite 
the tendency of the electric-vector linear polarization angle to 
align with the jet position angle found at radio frequencies on 
quasars. This implies that the sources in our sample do not 
contain the conditions required for their linear polarization an- 
gles to lie either parallel or perpendicular to the jet axis (in the 
observer's frame), i.e., that they have a markedly non axisym- 
metric character. Indeed, the 3 -dimensional character of a non- 



negligible populatio n of jets in AGN at different scales is known 



to be complex (e.g., Fomalont et al. 


2000c iPollack et al.l2003 


Jorstad et alj|2005l iLister & Homan 


l2005HA2udo et al.ll2006a 


2007). 



In the case of purely axisymmetric jets, cancellation of lin- 
early polarized emission parallel and perpendicular to the jet 
axis forces y to align only either parallel, or perpen dicular to 
the axis (e.g, iLvutikov et ai]|2005l; ICawfhorne)l2006l) . Our re- 
sults contradict this phenomenology, which suggests that the 
theoretical requirements for relativistic jets to show such behav- 
ior are not satisfied when observed at 3 mm. This implies that 
the 3 mm emitting region in relativistic jets in blazars is primar- 
ily non axially-symmetric at this wavelength. Hence, either the 
magnetic field or the emitting particle distributions (or both) 
responsible for the synchrotron radiation at 3 mm in a large 
fraction of blazars must have a markedly 3 dimensional (non- 
axisimmetric) character in order to account for the |y - 0j et | dis- 
tributions of our observations. 

We have shown that the circularly polarized emission de- 
tected from our observations is most likely generated by in- 
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trinsic synchrotron emission. This has an important conse- 
quence. Intrinsic circular polarization cannot be produced in 
pure pair (electron-positron, e + e~) plasma jets. Thus, sources 
in which such mechanisms govern circular polarization produc- 
tion must be composed, at least in part, by an electron-proton 
(p + e~) plasma. This may be taken as a tool for AGN jet compo- 
sition diagnostics through future mm wavelength polarimetric 
surveys. If a large fraction of sources show mc with moduli as 
large as those detected at cm wavelengths, this may imply that 
their mm jet emitting regions are primarily composed of p + e~ 
plasma. 
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DA 354 


13 57 04.4305 


+ 19 19 07.251 


0.719 


Q 


16.0 


146 


15 


9 




1406-076 




14 08 56.4811 


-07 52 26.665 


1.493 


Q 


19.6 


265 


43 


2 




1413+135 




14 15 58.8108 


+ 13 20 23.601 


0.247 


B 


20.5 


247 


15 


g 





12 
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TABLE 1 — Continued 



Source name 




— 

R.A. 


— 

Dec. 




Opt 


— — — 

V 


<»ia 
1° 1 


fobs 


Rel . 


LBAS 


(IAU) 


Alias 


(J2000.0) 


(J2000.0) 


z 


CI. 


mag 


[GHz] 


(4>ja) 


n - 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1418+546 




14 19 46.5784 


+54 23 14.616 


0.152 


B 


15.7 


133 


15 


9 




1502+106 




15 04 24.9752 


+10 29 39.080 


1.839 


Q 


18.6 


131 


86 


6 


Y 


1504-166 




15 07 04.7876 


-16 52 30.238 


0.876 


Q 


18.5 


205 


15 


8 




1510-089 




15 12 50.5321 


-09 05 59.845 


0.360 


Q 


16.5 


353 


86 


6 


Y 


1514-241 


APLib 


15 17 41.8190 


-24 22 19.431 


0.048 


B 


14.8 


173 


15 


9 


Y 


1546+027 




15 49 29.4326 


+02 37 01.069 


0.412 


Q 


17.4 


183 


86 


6 




1548+056 




15 50 35.2658 


+05 27 10.400 


1.422 


Q 


19.5 


1 


86 


6 




1606+106 




16 08 46.1974 


+10 29 07.666 


1.226 


Q 


18.7 


336 


15 


8 




1611+343 




16 13 41.0330 


+34 12 47.707 


1.401 


Q 


18.1 


164 


15 


8 




1622-297 




16 26 06.0237 


-29 51 26.770 


0.815 


Q 


19.5 


252 


43 


2 




1633+382 


4C 38.41 


16 35 15.4848 


+38 08 04.423 


1.807 


Q 


18.0 


288 


15 


8 


Y 


1637+574 




16 38 13.4457 


+57 20 23.874 


0.751 


Q 


16.9 


225 


86 


6 




1638+398 


NRA0 512 


16 40 29.6235 


+39 46 45.979 


1.666 


Q 


19.4 


282 


15 


8 




1641+399 


3C 345 


16 42 58.8001 


+39 48 36.958 


0.594 


Q 


16.6 


275 


15 


8 




1642+690 




16 42 07.8336 


+68 56 39.698 


0.751 


Q' 


20.5 


164 


15 


5 




1652+398 


Mrk 501 


16 53 52.2167 


+39 45 36.609 


0.033 


B 


13.8 


125 


15 


5 


Y 


1655+077 




16 58 09.0340 


+07 41 26.852 


0.621 


Q 


20.0 


333 


86 


6 




1657-261" 




17 00 53.1591 


-26 10 51.478 




u 












1716+686 




17 16 13.9209 


+68 36 38.684 


0.777 


Q 


18.5 


326 


5 


1 1 




1730-130 


NRAO 530 


17 33 02.7019 


-13 04 49.502 


0.902 


Q 


19.5 


31 


15 


8 




1732+389 




17 34 20.5664 


+38 57 51.398 


0.976 


Q 


20.6 


89 


5 


3 




1739+522 




17 40 36.9634 


+52 11 43.410 


1.379 


Q 


18.7 


41 


86 


6 




1741-038 




17 43 58.8510 


-03 50 04.604 


1.057 


Q 


20.4 


237 


86 


6 




1749+096 


OT081 


17 51 32.8104 


+09 39 00.700 


0.320 


B k 


16.8 


6 


15 


8 


Y 


1800+440 




18 01 32.2950 


+44 04 21.849 


0.663 


Q 


17.9 


271 


86 


6 




1803+784 




18 00 45.6222 


+78 28 04.022 


0.680 


B k 


15.9 


224 


86 


6 


Y 


1807+698 


3C 371 


18 06 50.6518 


+69 49 28.089 


0.050 


B 


14.2 


264 


15 


5 




1823+568" 


4C 56.27 


18 24 07.0480 


+56 51 01.484 


0.664 1 


B 


18.4° 


194 


86 


6 




1828+487 


3C 380 


18 29 31.8047 


+48 44 46.496 


0.692 


Q 


16.8 


316 


86 


6 




1830-211 




18 33 39.9093 


-21 03 40.049 


2.507 


Q 


18.7 










1842+681 




18 42 33.7085 


+68 09 25.034 


0.475 


Q 


18.1 


1 17 


86 


6 




1908-201 




19 11 09.6517 


-20 06 54.989 


1.119 


Q 


18.1 


4 


15 


9 


Y 


1921-293 




19 24 51.0545 


-29 14 29.838 


0.352 


Q 


18.2 


335 


86 


6 




1923+210" 




19 25 59.5932 


+21 06 26.106 




U 


16.1° 


245 


86 


6 




1928+738 




19 27 48.4595 


+73 58 01.592 


0.303 


Q 


16.1 


156 


15 


8 




1954+513 




19 55 42.7230 


+51 31 48.585 


1.223 


Q 


18.5 


306 


43 


7 




1957+405 


Cyg A 


19 59 28.3546 


+40 44 02.101 


0.056 


G 


15.1 


283 


15 


8 




1958-179 


20 00 57.0848 


-17 48 57.547 


0.652 


Q 


18.6 


207 


15 


8 




2005+403 




20 07 44.9340 


+40 29 48.622 


1.736 


Q 


19.0 


120 


15 


8 




2007+777 




20 05 30.9646 


+77 52 43.294 


0.342 


B 


16.7 


282 


86 


6 




2013+370 




20 15 28.7151 


+37 10 59.640 




B 1 " 


21.6° 


176 


86 


6 




2021+317" 


4C +31.56 


20 23 19.0066 


+31 53 02.395 




U 


19.0° 


168 


15 


8 




2023+336 




20 25 10.8256 


+33 43 00.265 


0.219 


B 




344 


86 


6 




2037+511 


3C418 


20 38 37.0188 


+51 19 12.687 


1.687 


Q 


21.0 


226 


86 


6 




2059+034 




21 01 38.8275 


+03 41 31.381 


1.015 


Q 


17.8 


30 


8 


10 




2113+293 




21 15 29.3850 


+29 33 38.540 


1.514 


Q 


20.6 


177 


15 


5 




2121+053 




21 23 44.4941 


+05 35 22.192 


1.941 


Q 


20.4 


276 


15 


8 




2128-123 




21 31 35.2540 


-12 07 04.725 


0.501 


Q 


16.1 


216 


15 


8 




2131-021 




21 34 10.3053 


-01 53 17.163 


1.284 


B k 


19.0 


100 


15 


8 




2134+004 


DA 553 


21 36 38.5791 


+00 41 54.319 


1.932 


Q 


17.1 


332 


15 


8 




2136+141 




21 39 01.3021 


+14 23 36.108 


2.427 


Q 


18.9 


310 


15 


8 




2145+067 


4C +06.69 


21 48 05.4509 


+06 57 38.710 


0.999 


Q 


16.5 


122 


15 


8 




2155-152 




21 58 06.2819 


-15 01 09.327 


0.672 


Q 


18.3 


207 


15 


8 




2200+420 


BL Lac 


22 02 43.2793 


+42 16 40.073 


0.069 


B 


14.7 


221 


86 


6 


Y 


2201+315 


4C 31.63 


22 03 14.9665 


+31 45 38.359 


0.298 


Q 


15.6 


209 


15 


8 




2210-257 




22 13 02.4963 


-25 29 30.054 


1.831 


Q 


19.0 


91 


5 


3 




2216-038 




22 18 52.0315 


-03 35 36.837 


0.901 


Q 


16.4 


188 


15 


8 




2223-052 


3C446 


22 25 47.2570 


-04 57 01.271 


1.404 


Q 


18.4 


128 


86 


6 




2230+114 


CTA 102 


22 32 36.4015 


+11 43 50.985 


1.037 


Q 


17.3 


126 


15 


8 


Y 


2234+282 




22 36 22.4627 


+28 28 57.525 


0.795 


Q 


19.0 


263 


15 


5 




2243-123 




22 46 18.2309 


-12 0651.110 


0.630 


Q 


16.4 


359 


15 


8 




2251+158 


3C 454.3 


22 53 57.7438 


+ 16 08 53.648 


0.859 


Q 


16.1 


255 


86 


6 


Y 


2254+617" 




22 56 17.9320 


+62 01 49.545 




u 












2255-282 




22 58 05.9656 


-27 58 21.312 


0.927 


Q 


16.8 


224 


15 


9 




2318+049 




23 20 44.8503 


+05 13 50.085 


0.623 


Q 


19.0 


322 


15 


9 




2345-167 




23 48 02.6085 


-16 31 12.022 


0.576 


Q 


18.4 


121 


86 


6 





Note. — Columns indicate, for each line: (1) IAU B1950.0 s ource name, (2) common sou rce name, (3) and (4) J2000.0 observ- 
ing right ascension and decl ination, respectively, (5) redshift from Veron-Cetty & Veron i 2006 ) catalog, (6) optical classification from 
Veron-Cetty & Veron (2006) catalog, (7) V magnitude from Veron-Cetty & Veron i 2006) catalog, (8) jet position angle (0j et ), (9) VLBI 
observing frequency of the image from which 0j el was measured, (10) references from which 0j C i was obtained, (11) source in the LBAS 
list of sources (Y), or not (...). 



References. — (l) I Agudo et alj UMfo\ ( 2) B U- Blazar Gr oup Web Page ; (3) Fomalont et al. (2000); (4) Jorstad et al. ffifflb ; (5) 
Kellermann eTaT] 120041) ; (6) |Lee e t al. i 2008); (7) Lister 1 2001 ): (8) Lister & Homan . 20Q3); (9) MOJAVE Web Page; (10) USNO Radio 
Reference Frame Image Database Web Page; (11) Xu et al. 1 19351). 

"Not present in Veron-Cetty & Veron 1 2006) catalog. 

Redshift from Sowards- Emmerd et al. i 2005). 

°y-mag= W-mag from NED. 

''Classified as a possible galaxy in Veron-Cetty & Veron (2006) catalog. 
'V-mag^optical-mag from NED. 

f From lAgudo et alj|2T)07l) ; IAcosta-Pulido et alj lin press!) , 
SRedshift from lWills ifeWrsl fT9751) . 
h Redshift fmmlNilssrm Rt al II200KH. 
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Table 2 



Summary of observing results. 



Source 




tint 




'"L 


X 


nic 


name 


Epoch 


M 


|Jv] 


!'■;■] 


r i 


[%] 


(1) 


(2) 


(3) 


(4) 


(5) 


16) 


(7) 


0003+380 


Jul. 2005 


4.0 


1.24±0.06 


2.20±0.54 


36.8± 7.0 


<0.67 


0048-097 


Jul. 2005 


8.0 


0.62±0.03 


< 1.72 




<0.89 


0059+581 


Sep. 2007 


8.0 


1.10±0.05 


< 1.58 




<0.73 


0106+013 


Jun. 2009 


4.0 


1.71±0.09 


< 1.58 




<0.63 


0112-017 


Jun. 2009 


16.0 


0.15±0.01 


<4.29 




<3.79 


0113-118 


Jul. 2005 


8.0 


0.96±0.05 


4.91±0.53 


157.0± 3.1 


<0.60 


0119+041 


Jul. 2005 


8.(1 


0.42±0.02 


< 2.00 




<1.41 


0133+476 


Jul. 2005 


4.0 


3.75±0.19 


4.11±0.51 


41.6± 3.6 


<0.44 


0135-247 


Jun. 2009 


12.0 


1.36±0.07 


3.70±0.55 


112.1±4.2 


<0.68 


0202+149 


Jul. 2005 


24.0 


0.42±0.02 


5.14±0.56 


114.1± 3.1 


<0.79 


0212+735 


Jul. 2005 


12.0 


1.51±0.08 


3.94±0.51 


144.0± 3.7 


<0.45 


0219+428 


Nov. 2008 


16.0 


().60±0.03 


4.36±0.54 


2.6± 3.6 


<0.73 


0221+067 


Jul. 2005 


32.0 


0.54±0.03 


6.62±0.54 


64.8± 2.3 


<0.70 


0224+671 


Jul. 2005 


8.0 


0.89±0.04 


4.29±0.54 


22. 1± 3.6 


<0.66 


0234+285 


Jul. 2005 


16.0 


3.87±0.19 


2.42±0.50 


126.7± 5.9 


<0.32 


0235+164 


Jul. 2005 


16.0 


1.96±0.10 


4.20±0.50 


14.7± 3.4 


<0.36 


0239+108 


Jul. 2005 


16.0 


0.25±0.01 


3.78±0.72 


27.8± 5.4 


<1.59 


0300+470 


Jul. 2005 


8.0 


0.45±0.02 


< 1.83 




<1.06 


0316+413 


Jul. 2005 


4.0 


5.95±0.30 


< 1.51 




0.46±0.11 


0333+321 


Sep. 2007 


8.0 


0.79±0.04 


1.99±0.57 


52.5± 8.3 


<0.79 


0336-019 


Jul. 2005 


24.0 


1.14±0.06 


4.20±0.51 


84.2± 3.5 


<0.43 


0355+508 


Jul. 2005 


4.0 


4.10±0.21 


2.05±0.50 


29.3± 7.0 


<0.34 


0403-132 


Jun. 2009 


16.0 


0.38±0.02 


4.01±0.63 


170.4± 4.4 


<1.25 


0415+379 


Jul. 2005 


4.0 


3.90±0.19 


< 1.52 




<0.37 


0420-014 


Jul. 2005 


12.0 


2.62±0.13 


2.70±0.50 


125.6± 5.3 


<0.35 


0422+004 


Jul. 2005 


12.0 


1.60±0.08 


5.52±0.51 


150.2± 2.7 


<0.43 


0430+052 


Jul. 2005 


20.0 


2.06±0.10 


< 1.51 




<0.35 


0439+360 


Jul. 2005 


4.0 


1.77±0.09 


< 1.55 




<0.55 


0454-234 


Jul. 2005 


8.0 


1.19±0.06 


4.51 ±0.54 


13.2± 3.5 


<0.71 


0458-020 


Jul. 2005 


24.0 


0.42±0.02 


2.78±0.58 


89.8± 5.9 


<0.87 


0514-161 


Feb. 2009 


8.0 


0.20±0.01 


<2.95 




<3.11 


0528+134 


Jul. 2005 


12.0 


2.00±0.10 


2.63±0.51 


91.2± 5.5 


<0.39 


0529+075 


Jul. 2009 


8.0 


0.98±0.05 


3.84±0.55 


149.2± 4.0 


<0.66 


0552+398 


May 2006 


8.0 


1.35±0.07 


4.72±0.54 


119.6± 3.3 


<0.64 


0605-085 


Jul. 2005 


8.0 


1.03±0.05 


4.50±0.54 


1.6± 3.4 


<0.70 


0607-157 


Jul. 2005 


8.0 


1.55±0.08 


3.12±0.52 


138. 8± 4.8 


<0.53 


0642+449 


May 2006 


8.0 


1.27±0.06 


2.75±0.54 


2.6± 5.7 


<0.66 


0716+714 


Aug. 2007 


8.0 


1.89±0.09 


11.44±0.53 


143.0± 1.3 


<0.65 


0727-115 


Jul. 2009 


8.0 


3.47±0.17 


2.55±0.50 


99.2± 5.7 


<0.35 


0735+178 


Feb. 2008 


8.0 


0.56±0.03 


2.56±0.61 


163.5± 6.8 


<1.00 


0736+017 


Oct. 2009 


12.0 


1.68±0.08 


< 1.52 




<0.38 


0745+241 


Jul. 2009 


8.0 


0.89±0.04 


5.72±0.54 


99.6± 2.7 


<0.62 


0754+100 


Jun. 2009 


16.0 


1.93±0.10 


5.44±0.51 


16.6± 2.7 


<0.64 


0804+499 


Jun. 2009 


12.0 


0.44±0.02 


2.24±0.60 


101.4± 7.7 


<1.01 


0805-077 


May 2006 


16.0 


0.89±0.04 


4.71±0.55 


175.2± 3.4 


<0.68 


0814+425 


Jul. 2009 


16.0 


1.03±0.05 


<4.76 




<2.74 


0820+560 


Feb. 2009 


8.0 


0.50±0.02 


2.94±0.63 


23.4± 6.1 


<1.09 


0823+033 


May 2006 


8.0 


1.92±0.10 


3.38±0.52 


27.5± 4.4 


0.60±0.16 


0827+243 


Feb. 2009 


8.0 


0.95±0.05 


< 1.60 




-0.63±0.20 


0829+046 


Aug. 2007 


16.0 


0.70±0.04 


8.20±0.62 


45.6± 1.9 


<1.21 


0834-201 


May 2006 


16.0 


1.02±0.05 


1.94±0.55 


43. 8± 8.0 


<0.79 


0836+710 


May 2006 


8.0 


2.53±0.13 


5.59±0.51 


107.6± 2.6 


-0.48±0.15 


0851+202 


Jul. 2005 


8.0 


3.20±0.16 


6.77±0.50 


152.8± 2.1 


<0.35 


0923+392 


May 2006 


8.0 


4.54±0.23 


4.16±0.50 


137.2± 3.5 


<0.34 


0945+408 


May 2006 


8.0 


0.60±0.03 


2.46±0.67 


47.0± 7.6 


<1.48 


0953+254 


Jul. 2005 


8.0 


0.90±0.05 


1.86±0.54 


41. 8± 8.3 


<0.65 


0954+658 


May 2006 


4.0 


2.33±0.12 


5.20±0.52 


8.1 ±2.9 


<0.82 


1012+232 


Jul. 2005 


8.0 


0.49±0.02 


3.85±0.61 


71.7±4.6 


<1.13 


1034-293 


Feb. 2009 


8.0 


1.20±0.06 


2.15±0.54 


138. 1± 7.3 


<0.61 


1039+811 


May 2006 


8.0 


0.45±0.02 


< 2.15 




<1.57 


1044+719 


Jul. 2005 


8.0 


1.12±0.06 


< 1.61 




<0.72 


1045-188 


Jul. 2005 


8.0 


1.05±0.05 


6.65±0.58 


148.1 ±2.4 


<0.86 


1055+018 


Jul. 2005 


8.0 


3.01±0.15 


4.01 ±0.50 


132.0± 3.6 


<0.35 


1116+128 


Feb. 2009 


8.0 


0.54±0.03 


5.02±0.63 


10.7± 3.7 


<1.22 


1124-186 


Jul. 2005 


8.0 


1.27±0.06 


9.56±0.53 


29. 5± 1.6 


0.58±0. 19 


1127-145 


Aug. 2007 


8.0 


1.31 ±0.07 


2.61±0.54 


81.9± 6.3 


<0.78 


1144+402 


Jul. 2009 


8.0 


0.77±0.04 


4.87±0.53 


117.2± 3.1 


<0.58 


1156+295 


Aug. 2007 


8.0 


0.72±0.04 


4.71±0.58 


161.4± 3.6 


<1.14 


1213-172 


Feb. 2009 


8.0 


0.93±0.05 


< 1.61 




<0.67 


1226+023 


Aug. 2007 


8.0 


16.01±0.80 


4.14±0.50 


131. 8± 3.5 


<0.31 


1228+126 


Jul. 2005 


4.0 


3.06±0.15 


7.82±0.52 


51.4± 1.9 


<0.55 


1244-255 


Jul. 2005 


8.0 


1.66±0.08 


2.50±0.52 


147.2± 5.9 


<0.52 


1253-055 


Aug. 2007 


8.0 


11.46±0.57 


< 1.50 




<0.31 


1308+326 


Jul. 2005 


24.0 


0.84±0.04 


6.18±0.54 


164.9± 2.5 


<0.69 


1328+307 


Oct. 2006 


24.0 


0.97±0.05 


13.40±0.52 


39.6± 1.1 


<0.45 


1334-127 


Jul. 2005 


4.0 


4.76±0.24 


7.13±0.50 


172. 1± 2.0 


-0.35±0.11 


1354+195 


Jul. 2005 


8.0 


1.54±0.08 


2.38±0.53 


21.8±6.4 


<0.58 


1406-076 


Aug. 2007 


8.0 


0.57±0.03 


< 1.93 




<1.62 


1413+135 


Sep. 2009 


16.0 


0.31 ±0.02 


< 1.84 




<U1 
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TABLE 2 — Continued 



Source 




tint 


Ss6 


m L 


X 


mc 


name 


Epoch 


[sj 


My] 


[%] 


[°] 


[%] 


a) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 A 1 O i IZA£. 

1418+546 


T..1 OAAC 

Jul. 2U05 


. 

1 2.0 


A OA_L_A A A 
0.80 ±0.04 


O OA_LA CA 

2.80±0.56 


1 A^7 H _L C O 

10/. /± 5. / 


<0.85 


i i i aa 
1 3UZ+ 1 wo 
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Note. — Columns indicate, for each line: (1) IAU B1950.0 source name, (2) observing epoch, (3) 
integration time, (4) 86 GHz flux density, (5) fractional linear polarization, (6) linear polarization electric 
vector position angle (7) fractional circular polarization. 
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FIG. 1 . — Sky distribution of the source sample in J2000.0 equatorial coordinates. 
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FIG. 2. — 86 GHz luminosity as a function of redshift. The dashed line indicates the luminosity for observer's frame flux density . 
line is for Sge = 5 Jy. Diamonds symbolize quasars, triangles denote BL Lacs, and squares correspond to radio galaxies. 
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FIG. 3. — Distribution of 86 GHz luminosity for all, quasars, BL Lac sourc es, radio galaxies , LBAS sources (i.e., those sources in our entire sample that were 
detected by Fermi-LAT during its first 3 months of operation (Tables 1 and 2 of Abdo et al. 2009)), and non-LBAS sources. Numbers in parentheses denote sample 
sizes of sources with known redshift. 
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FIG. 4. — Distribution of 86GHz luminosity for LBAS quasars, LBAS BL Lacs, non-LBAS quasars, and non-LBAS BL Lacs in our source sample. Numbers in 
parentheses denote sample sizes of sources with known redshift. 
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FIG. 5 . — Distribution of 15 GHz to 86 GHz spectral indices (»i5 86) for all sources in both the MOJAVE and our sample, and their corresponding quasar, BL Lac, 
LBAS, and non-LBAS sub-samples. The 15 GHz total flux density was taken from integrated intensity of MOJAVE images. For each source, the closest 15 GHz 
observation to our 86 GHz measurement was selected. Numbers in parentheses denote sample sizes. 
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FIG. 6. — Distribution of 86GHz fractional linear polarization for: all sources in the entire sample, quasars, high polarization quasars (HPQ), low polarization 
quasars (LPQ), BL Lac objects, radio galaxies, LBAS, and non-LBAS sources (from top to bottom). N is the number of sources in 1 % wide bins. Unfilled areas 
correspond to non-detections. Note that since we have chosen to consider 3<r upper limits with a,„ L ks 0.5 %, no data with ml ~ 1.6 % are available. Numbers in 
parentheses denote sizes of detected «?l samples. 
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FIG. 7. — Distribution of 86 GHz to 15 GHz fractional linear-polarization ratio for sources with detected linear polarization both in our survey and in MOJAVE. 
Two sources with m^/m^ 15 > 17 are not presented. The 15 GHz l inear polarization fractio n was computed from measurements of integrated total flux density and 
linearly polarized flux density from 15 GHz VLBI images listed in Lister & Homan 1 2005). Numbers in parentheses denote sample sizes. 
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FIG. 8. — Distribution of 86 GHz to 15 GHz fractional linear-polarization ratio for sources with detected linear polarization in our survey and th ose with detected 
core polarization in MOJAVE (i.e., m^/m^ 1 ^). Two sources with m^/m^^ > 14 are not presented. The used '«™J| measurements are those given by Lister & Homan 
1 2005). Numbers in parentheses denote sample sizes. 
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FIG. 9. — 86 GHz fractional linear polarization for sources with known redshift for the entire source sample, the quasar, the BL Lac, LBAS, and the non-LBAS 
sub-samples (from top to bottom). Arrows symbolize oil upper limits. The continuous line in the upper plot symbolizes the result of a linear regresion. Numbers in 
parentheses denote sample sizes. 
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FIG. 10. — Distribution of misalignment between the 86 GHz linear-polarization electric-vector position-angle (x, given in Tablef2) and the jet structural position 
angle (^>j e t, given in TablefT}. We present, from top to bottom, the entire source sample and the subsamples of quasars, BL Lacs, LBAS, and non-LBAS sources. 
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FIG. 12. — Distribution of \Q\ (top) and \U\ (bottom) Stokes parameters for sources with detected linear polarization. Four sources with \Q\ and/or \U\ in the 
range [300, 1 000] mJy are not shown. 
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FIG. 13. — Distribution of the absolute value of circular polarization for all major samples (labeled for each sub-plot) considered in this study. Black areas corre- 
spond to ttiQ detections at > 3<r. Grey shaded areas indicate observing results with > 2cr, whereas unshaded areas symbolize all niQ measurements, independently 
of their significance. 



